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INTRODUCTION 


The Fatio’s vole Microtus (Terricola) multiplex is probably pre- 
sent in Italy since last upper Pleistocene, despite the fact that its 
identification from other species of subgenus Terricola is not easy 
(CAPIZZI & SANTINI, 1999). In Italy, the present-day distribution is 
not totally known. M. ( T:) multiplex is considered to be an alpine and 
north-apennine species (LANZA & AZZAROLI, 1970). It lives scattered 
through the Alps, the north of the Apennines and as far as the north 
of Tuscany (SANTINI 1978; Capizzi & SANTINI, 1999) (Fig. 1). The 
species might be present also in the central Apennines (in Abruzzo) 
(AMORI et al., 1984). 


Generally M. (T.) multiplex is not found in mediterranean 
biotopes (AMORI et al., 1984), but some specimens were found in 
typical mediterranean stations of Lunigiana like Porto Venere (La 
Spezia) (SANTINI & Farina, 1978) and Pegazzano (La Spezia) (Rosi, 
2000). In these stations, M. ( T.) savii (de Sélys Longchamps) was not 
found, even though this species is present, on the contrary of M. ( T.) 
multiplex, at few kilometers in the lowlands near the cost (SANTINI & 
FARINA, 1978). 


This work is based on the study of different samples of M. ( T.) 
multiplex found in pellets of barn owl Tyto alba (Scop.) from lunigia- 
nese stations in provinces of Massa and La Spezia (Fig. 2). 


(*) UMR CNRS 5561 Biogéosciences et Laboratoire de Paléobiodiversité et de 
Préhistoire de l'Ecole Pratique des Hautes Etudes, Centre des Sciences de la 
Terre, Université de Bourgogne, 6 Bd. Gabriel, F-21000 Dijon (France) 

(**) Via Canarecchia 5, I-19021 Monti di Arcola (La Spezia, Italy). 
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The Ligurian stations are close to the sea cost, and they have a wet 
mediterranean climate: Pegazzano (SP) 150 m a.s.l.; Arcola (SP) 70 m 
a.s.l.; Isola (SP) 275 m a.s.l. The average annual temperature is 15.4°C 
in La Spezia, with an annual precipitation of 1404 mm (GAUDENZO 
et al., 1998). Way far from the sea, in the stations of Lunigiana, the 
temperature decreases and the precipitation progressively raises (Pon- 
tremoli, annual average temperature=12.4°C, humidity=71% and pre- 
cipitation=1667 mm, according to M. Ratti, pers. data); the stations 
of Lusuolo (MS) 179 m a.s.l., Malgrate (MS) 251 m a.s.l., Groppoli 
(MS) 363 m a.s.l., Pontremoli (MS) 236 m a.s.l. and Gravagna san 
Rocco (MS) 690 m a.s.l. are in this region. The stations of Fivizzano 
(MS) 326 m a.s.l. and Mazzola (MS) 449 m a.s.l. are in the Aulella 
valley, on the east part of Lunigiana. 


MATERIALS 


We examined first lower molars (M,) of 10 populations from Italy 
(denoted ITA): Pontremoli (Massa, Italy), n=58, Gravagna san Rocco 
(Massa, Italy), n=9, Pegazzano (La Spezia, Italy), n=14, Arcola (La 
Spezia, Italy), n=8, Isola (La Spezia, Italy), n=6, Malgrate (Massa, 
Italy), n=20, Lusuolo (Massa, Italy), n=14, Groppoli (Massa, Italy), 
n=6, Mazzola (Massa, Italy), n=48, Fivizzano (Massa, Italy), n=12. 


The materials have been collected in the years 1993-1996 and 
belong to the Rosi collection (Arcola, La Spezia, Italia); a part of them 
is deposited in the Centre des Sciences de la Terre, Université de 
Bourgogne (Dijon, France). 

The 10 populations were grouped according their respective geo- 
graphical distribution (Figure 2): group A: Pontremoli and Gravagna 
San Rocco; group B: Pegazzano, Arcola and Isola; group C: Malgrate, 
Lusuolo and Gròppoli; group D: Mazzola and Fivizzano. 

We compared Italian populations with: 

- 3 populations from Briangonnais (denoted BRI): Col du Lauta- 
ret (Hautes-Alpes, France) n=48, Col du Montgenèvre (Hautes-Alpes, 
France) n=17, Les Vigneaux (Hautes-Alpes, France) n=10; 

- 7 populations from Prealps of the Vercors, Trièves and around 
(denoted PRA): La Chapelle-en- Vercors (Drôme, France) n=12, Cha- 
beuil (Dróme, France) n=20, Marsaz (Dróme, France) n=3, Hauteri- 
ves (Dróme, France) n=2, Saint- Martin-de-la-Cluze (Isére, France) 
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Fig. 1 - Present-day distribution of Microtus (Terricola) multiplex in Italy. 
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Fig. 2 - Localities of samples of Microtus (Terricola) multiplex populations in Luni- 
giana (Italy): Pontremoli (Massa) (denoted 1), Gravagna San Rocco (Massa) (2), 
Pegazzano (La Spezia) (3), Arcola (La Spezia) (4), Isola (La Spezia) (5), Malgrate 
(Massa) (6), Lusuolo (Massa) (7), Gróppoli (Massa) (8), Mazzola (Massa) (9), 
Fivizzano (Massa) (10). 


n=10, Saint-Romans (Isère, France) n=24, Chantesse (Isère, France) 
n=14; 

- 7 population from Ticino and Wallis (denoted TIC): Gudo 
(Ticino, Switzerland) n=12, Bedano (Ticino, Switzerland) n=4, Isone 
(Ticino, Switzerland) n=4, Merde (Ticino, Switzerland) n=16, Biog- 
gio (Ticino, Switzerland) n=16, Varenzo (Ticino, Switzerland) n=11, 
Zermatt (Wallis, Switzerland) n=12. 

The populations from Ticino, assigned to the nominative subs- 
pecies multiplex (Fatio, 1905), described from Lugano (Ticino, Swit- 
zerland), and the population from Zermatt, assigned to the subspecies 
fatioi (Mottaz, 1909), described from Zermatt (Wallis, Switzerland), 
are pooled in the same geographical group due to the small sample size 
(n=12) of the population from Zermatt; 
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. The populations from Briançonnais were assigned to the sub- 
species druentius (Miller, 1912), described from Terres-Plaines near 
Barcelonnette (Alpes-de-Haute-Provence, France); 


The populations from Prealps and Triéves were assigned to 
the subspecies niethammeri (Brunet-Lecomte et Volobouev, 1994), 
described from Saint-Martin-de-la-Cluze (Isére, France), locality 
included in this geographical sample. 


METHODS 


Twenty three variables, denoted V1 to V23 (Figure 3), were 
taken on the occlusal surface of the M, (BRUNET-LECOMTE, 1988). 
We calculated 11 variables: the development of the anterior part 
AP=(V6-V3)/V6*100, the tilt of the pitymyan rhombus PR=V4-V3, 
the closure of the anterior loop AL=(V20-V18)/V21*100, the ratio 
length/width of M, V6V21=V6/V21 and the 7 other variables of the 
anterior part of the M, V10V9=(V10-V9)/V6*100, V11V9=(V11-V9)/ 
V6*100, V12V10=(V12-V10)/V6*100, V12V11=(V12-V11)/V6*100, 
V18V17=(V18-V17)/V21*100, V19V17=(V19-V17)/V21*100 and 
V20V19=(V20-V19)/V21*100. 
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Fig. 3 - Morphometry of the first lower molar of Microtus (Terricola) multiplex. 
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The comparative morphometry of M, was running on the total 
length of the M, TL=V6 and the 11 calculated variables AP, PR, AL, 
V6V21, V10V9, V11V9, V12V10, V12V11, V18V17, V19V17 and 
V20V19. The precision of measures was 10? mm, the variables TL 
and PR were expressed in mm, the other variables were without unit. 


The one-way analysis of variance (ANOVA) was used to compare 
the 12 criteria among Italian groups and materials from Briançonnais, 
Prealps-Triéves and Ticino-Wallis, completed by the Scheffé's test 
and the Dunnett's test if necessary. 


RESULTS AND DISCUSSION 

The morphology of the M, of Italian populations was visually 
similar to those of other populations of M. (T.) multiplex, characteri- 
sed by a well developed anterior part of the M, (AP>50) and a closed 
anterior loop (AL<30). 

Three sets of ANOVA were run (Table 1): 


1 — between the 4 geographical regions: 7 variables were signifi- 
cantly different: TL, AP, PR, AL, V6V21, V12V11 and V20V19. 


2 — between the 4 Italian groups: 6 variables were significantly dif- 
ferent: TL, V10V9, V11V9, V12V10, V12V11 and V20V19. 


Tab. 1 
ANOVA 


Region Group-t TIC 
«0.0001 i <0.0001 
<0.0001 0.2151 
<0.0001 (8777 
«0.0001 0.1471 
«0.0001 0.0113 

0.5424 alll 0.0046 

0.0181 Pa | 0.0008 


0.0955 0.0728 

«0.0001 0.0044 

0.4702 0.0864 

0.7472 0.7060 

VZONAO <0.0001 0.0362 


Tab. 1 - Analysis of variance (probability) for each of the 3 factors : geographical regions 
(denoted Region), Italian groups (denoted Group) and Italian groups of Italy com- 
pared to the Ticino and Wallis sample (denoted Group+TIC). 


COMPARATIVE MORPHOMETRY IN MICROTUS 465 


3 — between the 4 groups of Italy compared to Ticino and Wallis 
regions: 6 variables were significantly different: TL, V6V21, V10V9, 
V11V9, VI2VI1.and V 20V 129. 


ANOVA AMONG THE GEOGRAPHICAL REGIONS 


The Scheffé's test (Table 2) showed that the Italian sample was 
significantly different (p«0.05) from: 


- the Briangonnais sample for 6 criteria: TL (mean=2.787 
mm versus 2.433 mm), AP (mean=52.117 versus 51.743), V6V21 
(mean-2.632 versus 2.570), V11V9 (mean-1.120 versus 1.677), 
V12V11 (mean=16.336 versus 15.573) and V20V19 (mean=12.983 
versus 17.587); 


- the Prealps-Triéves sample for 5 criteria: TL (mean-2.787 mm 
versus 2.620 mm), AP (mean=52.117 versus 50.145), PR (mean=0.040 
mm versus —0.006 mm), V6V21 (mean=2.632 versus 2.560) and 
V12V11 (mean=16.336 versus 15.548); 


- the Ticino and Wallis sample for none criterion. 


Italian sample was characterised by a large size of the M, (mean 
TL=2.787 mm), a well developed anterior part (mean AP=52.117), a 
closed anterior loop (mean AL=20.332), a small difference for V11-V9 
(mean V11V9=1.120) and V20-V19 (mean V20V19=12.983), and a 
large difference for V12-V11 (mean V12V11=16.336). 


Summarising, the morphometry of the M, of Italian sample was 
generally similar to the Ticino and Wallis sample. 


ANOVA AMONG THE ITALIAN GROUPS 
The Scheffé's test (Table 3) showed that: 


- the group A was significantly different (p<0.05) from the group 
B for 2 criteria: V11V9 (mean=1.298 versus 0.476) and V20V19 
(mean=14.138 versus 10.717), and from the group D for 1 criterion: 
TL (mean=2.811 mm versus 2.720 mm); 


- the group B was significantly different (p<0.05) from the group 
A for 2 criteria: V11V9 and V20V19, from the group C for 1 criterion: 
TL (mean=2.895 mm versus 2.774 mm), and from the group D for 3 
criteria: TL (mean=2.895 mm versus 2.720 mm), V11V9 (mean=0.476 
versus 1.233) and V12V11 (mean=16.997 versus 15.852 mm); 
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"ab: 2 
TL (inmm) 


Region N Mean Std Dev Minimum Maximum ` Scheffé's test 


BENED E MM _ uo o 
TE 75 2.721 0.213 2.180 3.030 A 


BRI 75 2.433 0.189 2.070 2,190 e 


PRA 85 2.620 0.163 2.300 3.010 


B 
ITA 195 2.757 0.193 2.410 3.170 A 


AP (no unit) 


Region N Mean Std Dev Minimum Maximum ` Scheffé's test 
| —T7 ee MAIS 


ITA 195 52 INN 1.389 48.289 55.556 A 


PR (in mm) 


Mean Std Dev Minimum Maximum Scheffé’s test 


RIE 75 0.032 0.036 -0.060 0.100 A 
15 0.042 0.050 -0.090 0.150 A 
85 B 


PRA -0.006 0.036 -0.150 0.080 


195 0.040 0.037 -0.070 0.150 


AL (no unit) 


Region N Mean Std Dev Minimum Maximum ` Scheffé's test 


ene Mim 
TE 75 21.176 8.023 8.421 48.485 B 
BRI 75 25.852 9.871 10.227 58.696 A 


PRA 85 18.423 5.834 0.893 30.097 B 
B 


ITA 195 20:332 7.419 5.607 Ae AE 


V6V21 (no unit) 


Region N Mean Std Dev Minimum Maximum ` Scheffé's test 


85 : A A " B 
A 


I 
PRA 2.560 0.131 2-159 29950) 


ITA 195 2032 0.142 . 25262 3.326 
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V11V9 (no unit) 


Region N Mean Std Dev Minimum Maximum  Scheffé’s test 


uc 75 1309 [1.979 -1.394 4,545 AB 


BRI 75 677 10571 -1.327 53510. A 
PRA 85 1.361 1238 -3.175 4.603 AB 
ITA 195 1.120 {236 -2.062 4.377 B 


V12V11 (no unit) 


Region N Mean Std Dev Minimum Maximum Scheffé’s test 


TIE 75 Tosi 1.647 11.864 20.229 A 
BRI 75 15578 1747 11.024 19.485 B 
PRA 85 15.548 1597 12.346 20.238 B 


ITA 195 16.336 1.545 12.100 21.649 A 


V20V19 (no unit) 


Region N Mean Std Dev Minimum Maximum Scheffé’s test 
TIC 73 14.136 6.578 1.053 28.947 B 
BRI iS 17.587 4.959 8.333 27.835 A 

85 


PRA 11.060 3.948 0.000 21.698 C 
ITA 195 12.983 5:905 0.185 25.455 B.C 


Tab. 2 - Descriptive statistics and Scheffé's test of the criteria significant for the factor geo- 
graphical region. Scheffé's test: samples with the same letter are not significantly 
different (p>=0.05). 


- the group C was significantly different (p«0.05) from the group 
B for 1 criterion: TL, and from the group D for 2 criteria: V10V9 
(mean-14.709 versus 13.696) and V12V11 (mean-16.817 versus 
15.852); 

- the group D was significantly different (p«0.05) from the group 
A for 1 criterion: TL, from the group B for 3 criteria: TL, V11V9 and 
V12V11, and from the group C for 2 criteria: V10V9 and V12V11; 

- none difference was found between the groups A and C. 

The comparison with the sample from Ticino and Wallis as a refe- 
rence by the way of the Dunnett's test precised that: 

- the group A was significantly different (p«0.05) from the Ticino 
and Wallis sample for 1 criterion (TL: mean = 2.811 mm versus 2.720 
mm); 
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Tab. 3 


TL Gn mm) 


2:210 0.149 2.410 3.030 


V6V21 (no unit) 


Group N Mean Std Dev Minimum Maximum Dunnett’s test 


EE CU E A. ue. __ 
mean TIC-2.582 
A 67 2.633 0:132 2.420 2.932 NS 


21629 0.134 2:920 3.000 NS 


V10V9 (no unit) 


bh W —___r__e__r—_—»m lonsòi@_@ i cane RR 
pee emmmer mg o e... e —_onnniir 
A 67 14.616 1.785 9.854 17.871 AB NS 
B 28 14.439 0.935 12.542 16.140 ‘AB NS 
C 40 14.709 1.603 11.661 185315 A NS 
60 B S 


D 132696 1:695 8.897 17.143 


Group . N Mean Std Dev Minimum Maximum Scheffé’s test Dunnett’s test 


V11V9 (no unit) 


LAOS 1:289 -2.062 3.0198 


COMPARATIVE MORPHOMETRY IN MICROTUS 469 


V12V10 (no unit) 


-0.755 6.709 


V12V11 (no unit) 


67 16.206 1.454 13.028 19.870 


28 163087 1.653 13.448 201139 A 
40 16:817 1.347 14.583 20.619 A 


157852 15553 12.106 21.649 B 


V20V19 (no unit) 


iena NUM ___ —m_e | 
A 67 i . 58d 1.852 — 24742 A NS 
B 28 —— 10.717 4.397 3.696 20.561 B S 
E 40 12.236 4.593 3.571 25.000 AB NS 
D 60 13.249 5.783 1.942 25455 AB NS 


Tab. 3 - Descriptive statistics, Scheffé’s test and Dunnett’s test of the criteria signifi- 
cant for the factor Italian group. Scheffé’s test: samples with the same letter 
are not significantly different (p>=0.05). 

Dunnett’s test: NS = sample not significantly different (p>=0.05) from 
Ticino and Wallis sample, S = sample significantly different (p<0.05) from 
Ticino and Wallis sample. 


- the group B was significantly different (p<0.05) from the Ticino 
and Wallis sample for 4 criteria (TL: mean = 2.895 mm versus 2.720 
mm, V6V21 : mean=2.682 versus 2.582, V11V9 : mean=0.476 versus 
1.369 and V20V19 : mean=10.717 versus 14.136); 

- the group C was significantly different (p<0.05) from the Ticino 


and Wallis sample for 1 criterion (V11V9: mean = 0.756 versus 
1.369); 
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- the group D was significantly different (p<0.05) from the Ticino 
and Wallis sample for 1 criterion (V10V9: mean = 13.696 versus 
14.574). 


Summarising, the group B was the most different from the Ticino 
and Wallis sample, and within the Italian samples, the groups B and D 
were the most separated, while the groups A and C occupied an inter- 
mediate position between the groups B and D. 


SYSTEMATIC POSITION OF THE ITALIAN POPULATION 


Previous studies (GRAF & MEYLAN, 1980, BRUNET-LECOMTE & 
VOLOBOUEV, 1994, HARING et al., 2000) have pointed-out the gene- 
tic and karyological polymorphism of the species. GRAF & MEYLAN 
(1980) studied by electrophoretic and cytogenetic analyses the popu- 
lations of Varenzo (Ticino) and Fivizzano (Massa), 2 samples belon- 
ging to the Ticino and Wallis sample (Varenzo) and to the sample ` 
D (Fivizzano) of the Italian group, reporting a very small genetic 
distance (Nei's distance = 0. 01) between the samples of Varenzo and 
Fivizzano; the sample from Fivizzano differed from the sample from 
Varenzo by the loss of a pair of acrocentric autosomes, while the sex 
chromosomes were similar to those from Varenzo sample. 


The morphometric analysis of the M, confirmed the close 
relationship between the group D (Fivizzano and Mazzola) and the 
Ticino and Wallis sample and showed a medium difference, in the size 
like in the shape, between the group D (Fivizzano and Mazzola) and 
the group B (the 3 samples of La Spezia). This difference is intere- 
sting to report, because of the geographical position of the group B 
near the sea. 


Recently HARING et al. (2000) have confirmed that the Tuscan 
populations belong to M. ( T.) multiplex. 


The synthesis of genetic, cytogenetic and morphometric 
data leads to propose the following hypothesis: because of the small 
genetic distance between the samples of Varenzo and Fivizzano, the 
karyotypic and the morphometric differences found between the Tici- 
nese and the Italian populations are recent. 


The available data lead to classify the Italian population in the 
nominative subspecies multiplex. 
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CONCLUSIONS 


The morphometric analysis of the M, confirmed the close rela- 
tionships displayed by genetic and cytogenetic analyses between Tici- 
nese and Lunigianese populations of M. ( T.) multiplex. On the other 
hand, this analysis showed a medium difference between the sample 
from the coast of the province of La Spezia (Pegazzano, Arcola and 
Isola) and the samples from the province of Massa, particularly the 
sample D (Mazzola and Fivizzano). At the moment it is not possible to 
know if this difference is correlated to genetic difference between the 
populations or has an environmental explanation. Next studies, like 
genetic analyses, will try to explain this difference. 
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ABSTRACT 


'The comparative morphometry of the first lower molar of different populations 
of Microtus ( Terricola) multiplex from Lunigiana (Italy) confirmed that these popula- 
tions are close to the nominative subspecies multiplex from 'Ticino (Switzerland) and 
showed that, within the Lunigianese populations, the population of the province of la 
Spezia is the most differentiated from the nominative subspecies. 


RIASSUNTO 


La morfometria comparativa del primo molare inferiore delle diverse popolazioni 
di Microtus (Terricola) multiplex (Fatio) della Lunigiana (Italia) ha confermato che 
queste popolazioni sono in stretta relazione con la sottospecie topotipica del Canton 
Ticino (Svizzera); inoltre, fra le popolazioni lunigianesi, quella della provincia della 
Spezia appare la piu differenziata dalla sottospecie nominale. 


